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Abstract 

The current weak lensing measurements of the large scale structure are mostly 
related to statistical study of background galaxy ellipticities. We consider a pos- 
sibility to extend lensing studies with intrinsically unresolved sources and suggest 
that spatial fluctuations in the integrated diffuse emission from these sources can 
be used for a lensing reconstruction. Examples of upcoming possibilities include 
the diffuse background generated by dusty starburst galaxies at far-infrared wave- 
lengths, first stars and galaxies in near-infrared wavelengths, and the background 
related to 21 cm emission by neutral gas in the general intergalactic medium prior 
to reionization. While methods developed to extract lensing information from cos- 
mic microwave background (CMB) temperature and polarization data can be easily 
modified to study lensing properties using diffuse backgrounds at other wavelengths, 
we suggest that the lensing extraction from these backgrounds using higher order 
non-Gaussian clustering information alone may not be the best approach. In con- 
trast to CMB anisotropics, reasons for this include the lack of features in the clus- 
tering power spectrum such that the resulting lensing modification to the angular 
power spectrum of low-redshift diffuse backgrounds, at arcminute angular scales, 
is insignificant. While the use of low redshift backgrounds for lensing studies will 
be challenging, due to confusing foregrounds among other reasons, the use of sug- 
gested backgrounds will extend the reconstruction of the integrated matter power 
spectrum out to redshifts of 15 to 30, and will bridge the gap between current and 
upcoming galaxy lensing studies out to, at most, a redshift of a few and planned 
weak lensing studies with CMB out to the last scattering surface at a redshift of 
1100. 



1 Introduction 



The weak gravitational lensing of background galaxies by the intervening large 
scale structure is now a well known probe of physical cosmology (Kaiser 1992; 
Kaiser 1998; Jain & Seljak 1997; Hu & Tegmark 1999; Bartelmann & Schneider 
2001). The current weak lensing studies are mostly limited to reduced-shear 
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estimates that involve statistical study of galaxy cllipticities or shape infor- 
mation reconstructed through expansions related to certain orthogonal basis 
functions (for a recent review, see, Refregier 2003). Such a lensing reconstruc- 
tion, unfortunately, is limited to resolved background sources for which reh- 
able shape measurements, after accounting for the point spread function, are 
possible. While statistics of current and most upcoming weak lensing surveys 
are expected to be limited by instrumental constraints, such as issues related 
to the point spread function, with reasonable improvements on the observa- 
tional side, the statistics related to large scale structure lensing studies will 
eventually be limited by the number of resolved background sources. In this 
limit, to further improve mass reconstruction related to lensing, background 
objects which are intrinsically point-sources, and, thus, for which shape mea- 
surements are impossible, must be included in the analysis. A situation where 
point-sources dominate resolved source number counts is already present in 
some of the deepest space-based images available to date, including North and 
South Hubble Deep Fields and recent images from the Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope. 

As one moves to higher source redshifts, beyond what is currently encountered 
in large scale structure lensing studies, the importance of point sources is more 
evident; The first stars and proto galaxies, at redshifts of order 15 to 30 and 
higher, provide an ideal background for lensing studies given the large path 
length to these background sources. While the detection of an individual first 
star is beyond the capabilities of even next generation instruments, the spatial 
fluctuations in the integrated diffuse background emission from these stars, at 
least in most optimistic models of their emission, may be within reach (e.g., 
Cooray ct al. 2003). During the transit to us, this diffuse emission is expected 
to be affected by the intervening mass distribution via gravitational lensing 
and provides a background for a potentially interesting lensing study. 

Contrary to background intensities at optical and infra-red wavelengths (Cooray 
et al. 2003; Knox et al. 2001), where the emission can eventually be broken 
down to point sources, one also encounter backgrounds which arc composed 
of truly diffuse emission that cannot easily be separated to individual objects. 
Examples of such diffuse backgrounds are the 21 cm rest wavelength emission 
from neutral Hydrogen in the intergalactic medium prior to reionization (Scott 
& Rees 1990; Madau et al. 1997; Tozzi et al. 1997; Iliev et al. 2002; Furlanetto 
et al. 2003) and, of course, the cosmic microwave background (CMB) radia- 
tion. The lensing effect on CMB temperature and polarization ffuctuations, 
that were generated at the last scattering surface (Peebles & Yu 1970; Sunyaev 
& Zel'dovich 1970; Silk 1968; Hu & Dodelson 2001), is now weU understood 
(Scljak 1996; Metcalf & Silk 1997; Hu 2000) and various techniques have been 
developed to extract lensing information from high resolution CMB tempera- 
ture and polarization maps (Zaldarriaga & Seljak 1999; Benabed et al. 2001; 
Guzik et al. 2000; Hu 2001; Hu & Okamoto 2002; Cooray & Kesden 2003; 
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Hirata & Seljak 2003). These techniques make use of the fact that at the last 
scattering surface, where primary CMB temperature and polarization fluctu- 
ations are generated, these anisotropies follow a Gaussian distribution, while 
with gravitational lensing deflections when propagating to us, the fluctuations 
pattern on the sky is non-Gaussian (Bernardeau 1997; Hu 2001; Zaldarriaga 
2000; kesden et al. 2002). The lensing extraction techniques make use of mea- 
surements related to certain non-Gaussian aspects, such as a collapsed form of 
the trispectrum (Zaldarriaga 2000; Hu 2001), which is related to a filtered ver- 
sion of the squared-temperature power spectrum (Cooray & Kesden 2003), to 
directly estimate and reconstruct the deflection angle or the projected lensing 
potential power spectrum. 

In the case of CMB, the study of lensing and its extraction is usually consid- 
ered under a perturbative approach with the deflection angle as a small change 
to the anisotropy on the sky towards a given direction (Hu 2000; Goldberg 
& Spergel 1998; Cooray & Hu 1998; Zaldarriaga 2000). Under such an ap- 
proach, for most practical purposes, the lensing effect can be calculated up 
to the second order in the deflection angle or first order in the lensing poten- 
tial power spectrum. Since the rms deflection angle for a photon propagating 
from a redshift of 1100 is of order a few arcminute, a resulting complication 
here is that perturbative description of the lensing effect breaks down at small 
angular scales at and below the rms deflection angle. This, however, is not 
a major concern for the CMB lensing description since there is no intrinsic 
small angular scale power at these same angular scales due to the presence of 
signiflcant damping. On the other hand, when compared to CMB tempera- 
ture fluctuations, diffuse backgrounds at lower redshifts are expected to show 
substantial clustering power at arcminute angular scales and the perturbative 
approach just to the first order in the deflection power spectrum may not be an 
adequate description. Another important difference, when compared to CMB 
anisotropies that were generated at the last scattering surface with narrow 
range in redshift, is that background sources may be distributed over a wide 
range in redshift and this distribution must be properly accounted as well. 

While there are complications from foregrounds and other confusions, we sug- 
gest that techniques developed to extract lensing information from CMB data 
can also be extended and applied at other wavelengths for lensing studies. 
In the case of CMB anisotropies, there are substantial modifications to the 
anisotropy angular power spectrum from the lensing effect due to the acoustic 
peak structure and the presence of a damping tail such that lensing transfers 
power from large scales to arcminute scales where the lensing effect is easily 
detectable. When compared to a lensing study with CMB, the lensing infor- 
mation one can extract from low redshift diffuse backgrounds is significantly 
limited. As we will discuss, this results from the fact that, at arcminute an- 
gular scales, the lensing modification to the power spectrum is minor and the 
non-Gaussianities generated are relatively smaller. This, again, is due to the 
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lack of features in the clustering power spectrum such that the modification 
to the clustering pattern, at arcminutc angular scales of interest for future 
observations, remains the same under gravitational lensing. The lack of signif- 
icant differences between the intrinsic and lensed power spectra, thus, limits 
the use of statistics such as clustering variance for lensing purposes. 

To extract lensing information better, one require arcsecond scale clustering 
information, though at such small angular scales, a substantial number of 
sources that contribute to the diffuse background will begin to be resolved. 
In the limit that most or all sources are resolved to point-like objects, instead 
of shear reconstruction based on source shapes, statistics related to the dis- 
tribution of the point sources can also be used for a shear estimation. As an 
example, instead of extracting information from the non- Gaussian clustering 
information, one can make use of the anisotropy of clustering generated by 
foreground lensing. This anisotropy can be used as an estimator of shear and 
can be approached in the same manner ellipticity information from galaxy 
shapes are used for an estimate of shear. Such a study has already been at- 
tempted towards known regions of mass concentrations, such as galaxy clusters 
(van Waerbeke et al. 1997). Similar studies can be extended to large angular 
scales and the large scale structure as a whole when carefully selected back- 
grounds are selected for lensing reconstruction so as to minimize confusions 
resulting from foregrounds. Another aspect is the change resulting from lens- 
ing magnification and several suggestions have already been made to include 
unresolved sources in weak lensing studies through variations associated with 
number counts and, when measurable, source sizes (Jain 2002; Pen 2003). 

The paper is organized as follows. In the next section, we discuss the lensing 
effect on diffuse backgrounds and extend calculations related to weak lensing 
effect on CMB anisotropies. As an example, we discuss the lensing modification 
of the near-IR anisotropy fiuctuations due to first stars following calculations, 
related to spatial clustering discussed in Cooray et al. (2003). We discuss the 
extent to which lensing information can be extracted from this and other back- 
grounds and a comparison to the lensing extraction with CMB data. We also 
consider potential biases and systematic effects for a lensing reconstruction 
and conclude with a summary. 



2 Calculation 

2.1 Lensing Effect on Clustering of Diffuse Backgrounds 

In order to derive the weak lensing effect on angular clustering properties 
of any diffuse background, we follow a calculation similar to that introduced 
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for the lensing extraction from CMB temperature anisotropies (e.g., Hu 2001; 
Seljak & Zaldarriaga 1999; Cooray & Kesden 2003) and make use of the flat sky 
approximation. When compared to CMB, there is one important difference in 
that one should properly account for the broad redshift distribution of sources 
that contribute to the background light, instead of the assumption related to a 
narrow source redshift distribution. Thus, we break the background intensity 
as 

oo 

/tot(n) = j drdl{rh,r) , (1) 





where 51{rh, r) is the fractional contribution to the total emission as a function 
of the radial distance or look-back time, from the observer: 

^ ' 



when the expansion rate for adiabatic CDM cosmological models with a cos- 
mological constant is given by 



Hi 



Q^(1 + z)^ + Qk(1 + -z)' + ^^a 



(3) 



Here, Hq can be written as the inverse Hubble distance today cHq^ — 2997.9/i~^ 
Mpc. We follow the conventions that in units of the critical density 3Hq/87iG, 
the contribution of each component is denoted = c for the CDM, b for the 
baryons, A for the cosmological constant. We also dcflnc the auxiliary quan- 
tities Qrn = flc + and fix = i — J2i^i, which represent the matter density 
and the contribution of spatial curvature to the expansion rate respectively. 
Note that, though we discuss a general derivation of lensing effect on back- 
ground source clustering, we show results for the currently favorable ACDM 
cosmology with flj, = 0.05, flrn = 0.35, fl\ = 0.65 and h = 0.65. 

Weak lensing deflects the path of background photons resulting in a remap- 
ping of the observed anisotropy pattern on the sky, such that the fractional 
contribution from each redshift is modified as 

(5/(rn, r) — 5I[r{h + V0(rn, r)}, r] 
^ SI{rh, r) Vi0(rn, r) V'(5/(rn, r) ^Vj0(rn, r) Vj0(m, r) V'V^'(5/(rn, r) . .(.4) 

Here, 5/(rn, r) is the unlensed fractional component of the background light 
from a distance of r, (5/(rn, r) is the same contribution when affected by the 
gravitational lensing deflections during the transit, and V0(rn, r) represents 
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the lensing deflection angle for a photon propagating from a distance of r. The 
second line of Eq. 4 treats the lensing deflection as a perturbative parameter. 
Typical lensing calculations related to CMB (e.g., Hu 2000), treat lensing ef- 
fect on the angular power spectrum to the second order in this perturbation 
expansion, while lensing extraction techniques related to CMB consider in- 
formation from the first order correction (e.g., Seljak & Zaldarriaga 2000; Hu 
2001). Such an approach works well since CMB has no intrinsic anisotropy 
at small angular scales below few arcminutes, though the typical deflection 
angle of a photon from the last scattering surface is ~ 2.5 arcminutes. The 
accuracy to which lensing can be extracted from CMB under this perturbative 
approach is discussed in Hirata & Seljak (2003). Here, for low redshift diffuse 
backgrounds, the perturbation expansion up to second order in the deflection 
angle is not adequate to fully account for the lensing effect on the cluster- 
ing power spectrum. We extend the calculation explicitly to higher order and 
present an exact derivation of the lensed power spectrum which should be 
valid even at small angular scales where the lensing deflection angle is of the 
same order. To properly evaluate the lensed contribution, however, one must 
evaluate an infinite series of integrals which is numerically exhaustive; here, in 
presenting numerical results, we only consider the lensed power spectrum up to 
the fourth order in deflection angle or the second order in the deflection-angle 
power spectrum. 

In equation 4, where V0 is the deflection angle, is a radial projection of the 
gravitational potential, $ (see, e.g. Kaiser 1992), for a source at a distance Vg 



Note that as Q,k — > 0, — > r. 

Taking the Fourier transform, as appropriate for a flat-sky, we write the mod- 
ification to the fraction contribution arising at a distance r as 





where the comoving angular diameter distance is 



(6) 




(7) 



where 
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L{h, W) ^ - 1/) (1, - 1/) ■ 1/ - ^ / 

X 0(1, - 1,' - 1/0 (1/' . i,')(ii - 1/ - h") -^^'+11^.1 

X 0(ii - 1/ - h" - W") ih" . . hO (U - ii' - ii" - ii'") • 1/ . + .... 

(8) 



Here, for simplicity, we have dropped the cxphcit dependence on the radial 
distance, r, in SI and 0, and have expanded the correction to the third order 
in (j). 

We define the power spectrum of the total intensity fiuctuation field, in the 
fiat sky approximation, following the usual way 



(7(1)7(10) ^(27rf5D(l + l')a, (9) 

where is the Dirac delta function. Similarly, the power spectrum of the 
lensed fluctuation fleld, with 7(1) replaced by 7(1), is Ci. 

Note that in the absence of lensing, the observed anisotropy power spectrum is 
simply the sum of anisotropy contributions over the normalized source redshift 
distribution such that Ci = J drCf^ . The observed angular power spectrum, 
however, consists of both the unlensed intensity and a perturbative correction 
related to the lensing effect. To calculate the final lensed-clustering power 
spectrum, we first substitute Eq. 7 in Eq. 1 and reintroduce the dependence 
on distance. The lensing contribution to the power spectrum is calculated 
following approaches used to describe lensing effect on CMB anisotropics (e.g., 
Hu 2000) and after some straight forward calculations, we obtain, to the fourth 
order in or second order in Cf: 



J drR{r)Cf^ 



1 - 2 



Ci^Ci- t- 

(10) 



+ 



where Cf^ is the fractional contribution to the intensity fluctuation power 
spectrum as a function of the radial distance (see. Section 2.2) and 



R{r) ^^j dll'Cfir) . 



(11) 
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The cumulative variance of intensity fluctuations is conserved under lensing 
such that (7^ = / (i^I/ (27r)^Q = / dH/ (27r)^C/; this is equivalent to the fact that 
lensing does not create or destroy power but rather results in a redistribution 
such that surface brightness remains conserved. Under lensing, however, the 
filtered variance of fiuctuations as a function of angular scale is different such 
that a'^{e) = JdH/{2nyCiW^{W), where W{W) is the Fourier transform of 
the filter — such as top-hat, Gaussian etc. — is different from that expected 
from the unlensed clustering power spectrum, / dH/ {27TyCiW^{l9). While the 
modification to variance resulting from lensing can be used for a reconstruction 
of convergence (Pen 2003), such a measurement only works out accurately 
if the differences between Ci and Ci is significant, as a function of angular 
scale or multipole /, under gravitational lensing. While this is the case for 
arcminute scale power spectrum of CMB anisotropics, we will show that the 
difference between intrinsic- and lensed-power spectra of diffuse backgrounds 
is insignificant and will complicate a straight forward lensing reconstruction 
based on the variance alone. 

In the case of CMB anisotropics, in Eq. 1, SI{rh, r) = P^^{n)5D{r — rg) where 
To is the distance to the surface of last scattering and P^^{n) is the total 
CMB intensity such that for spatial fiuctuations, Cf^ — Cf^^Soir — ro) with 
the primordial anisotropy power spectrum given by Ci^^. This assumption 
simplifies Eq. 10 to a result valid for any background source of narrow width 
in distance or in redshift space, as 



dH^ 



CMB 



l-/2i?(ro) + ^i?2(ro) 



^II^J 1^^^ (-o)c,t(-o)qiT-..|[(i - ii - h) ■ im - h - h) ■ h? . 



(2 



(12) 



This equation is equivalent to that in Hu (2000) for weak lensing effect on 
CMB anisotropy power spectrum when one drops higher order terms which 
are proportional to [Cf]^, and CfR. As we will discuss, to describe the 
lensing effect related to CMB, these higher order terms are not crucial. 

When one extends the perturbative calculation further, the power spectrum 
can be represented as a sum of a series of integrals such that 



^CMB _ ^CMB y^J^_-^y ^'^^Rj^oY 
1=0 ^' 
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(S^^^ (-o)...Qt(ro)q-^„,[(l - ... - IJ . - ... - IJ • 



(27r)2 7 (2 

|l-...-_l_„|2*i?(ro) 



xE(-i r ^ ' ^ (13) 

j=0 



where the last term is a generahzation of this series to the n-th term. The 
summations over R{ro) in each of these terms can be simphfied analytically 
and we obtain 



^CMB ^ ^CMBg-;2R(ro) 



(27r)2 
12 



2 



+ \lWf'l ^C-t('-o)C,t(ro)qiT-.,l(l - 1. - 1.) • - 1. - I2) ■ hf 

xe-|i-ii-i2l'«(^o) ^ ^ 

^hJ Wf- 1 |!^'^'t('-°)-c.t('-.)cf-"^,i(i - - - 1.) • 1.1^-1(1 - ... - 1„) • 1, 

xe-l'--'"l"«l'-"> . (14) 

Though this expression is useful since it presents the exact calculation related 
to how gravitational lensing modifies the clustering power spectrum, when 
numerically evaluating the lensed power spectrum with a limited number of 
terms in this series, one cannot simply evaluate the first few terms and ignore 
the rest. To evaluate the contribution accurately, one should include terms 
which are same order in Cf and R] For example, to the second order in Cf, 
contributions include three terms from the first line in Eq. 14, two terms from 
the second and a term from the third, when exponentials in each of the above 
lines are expanded. 

Note that the lensing effect only modifies the clustering pattern of an intrin- 
sically clustered background. If the spatial distribution of sources responsible 
for the diffuse background is simply Poisson, then one expects a shot-noise 
like angular power spectrum for fluctuations in the intensity such that Cf^ is 
a constant. In this case, Eq. 10 can be simplified and leads to Ci = Ci such 
that spatial fluctuations, after lensing, remains a Poisson distribution with the 
same shot-noise. This is again consistent with the fact that lensing does not 
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create or destroy power but only results in a modification to how power is 
distributed. 



The same can also be deduced from Eq. 14 by taking into account the exact 
calculation instead of the approximate one in Eq. 10. In the case of an intrinsic 

Poisson-clustered power spectrum, or a smooth or slowly varying intrinsic 
power spectrum, one can set C'|^^™_i^| = Q^mb ^j^g substitution that 
|1 — ... — 1„| ~ /. In this case, one can write Eq. 14 as 



^CMB _ ^CMBg-i2i?(ro)|-j^ _|_ 



(16) 



which simplifies as Cp^^ — Cp^^ since the series within curly brackets is 
simply e' This is important since as we discuss later, lensing only modifies 
the clustering pattern of background sources if the distribution has distinct 
signatures, such that the angular power spectrum show features. In the case 
of a smooth distribution, which may also be Poisson, lensing does not modify 
the clustering pattern significantly. In this sense, most shot-noise dominated 
backgrounds such as radio sources and galaxy clusters through the SZ effect 
at arcminute scale CMB experiments, will not be modified by gravitational 
lensing. We find a similar situation for some what higher redshift backgrounds 
at IR and other wavelengths. The only lensing modification is then associated 
with magnification changes to source number counts (Tegmark & Villumsen 
1997) such that, down to a certain fiux limit, more sources are magnified 
towards regions of high mass concentrations and the extent to which this is 
important depends strongly on the distribution of number counts as a function 
of fiux. For truly diffuse backgrounds such as CMB and 21cm emission from 
the general IGM prior to reionization, such an effect is not expected. 

In Eq. 10 and after, Cf (r^) is the power spectrum of lensing potentials, for 
a source at a distance of r^. In terms of the three-dimensional matter power 
spectrum, P{k), we can write this power spectrum of lensing potentials as 



Is Is 



where 
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Fig. 1. Angular power spectra of the lensing deflections as a function of the source 
redshift Zg. The sohd hnes show the fuhy non-hnear contribution to the power spec- 
tra, while the dotted lines show the contribution under linear theory power spec- 
trum. While low redshift lensing requires a fully non-linear treatment, for sources at 
redshifts of 15 and higher, non-linearities do not affect lensing deflections at angular 
scales at least of few arcminutes and below (Z ~ 10^). 

In terms of the well known convergence, power spectrum, that define the 
integrated surface density of mass out to the background distance r^, one can 
write Cf = I'^/ACf. In calculating the matter power spectrum, we use the 
fitting formulae of Eisenstein & Hu (1999) to evaluate the transfer function 
for CDM models and set the normafization to match both COBE and clusters 
(Bunn & White 1999; Viana & Liddle 1999). Note that in finear theory, the 
power spectrum can be scaled in time, P{k,r) = G^{r)P{k,0), using the 
the growth function (Peebles 1980). In the non-linear regime, one can use 
prescriptions such as the fitting function by Peacock & Doods (1996) or halo- 
based approaches (Cooray & Sheth 2002) to calculate the fully non-hnear 
density field power spectrum. 

We illustrate these power spectra as a function of Zg, the source redshift, in 
Fig. 1. Note that for sources at low redshifts, ~ 1, non- linear corrections 
to the matter power spectrum is important since lensing signal from degree 
angular scales and below are dominated by non-linearities. As one moves to 
a higher source redshift, non-linear corrections move to smaller angular scales 
such that for sources at redshifts 15 and higher, non-linear corrections are no 
longer important as lensing statistics at arcminute scales are fully described by 
linear fiuctuations in the intervening matter density field. For CMB photons. 
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the rms deflection angle, 6'rms = [/ (i^l/(27r)^/^Cf]°'^, is 2.6' while for a source 
at a redshift of 15, ^^ims ~ 2'. The fact that the typical deflection angle is 
of order few arcminutes, require that clustering modifications not be studied 
just to the first order perturbative correction and that higher order terms be 
included. 

In Fig. 2, we show the CMB anisotropy power spectrum and the resulting 
modification related to the gravitational lensing effect. Note that the lensing 
effect moves power from large angular scales, where acoustic peaks are found, 
to small angular scales (Seljak 1996; Hu 2000). This results in an increase of 
power at arcminute angular scales where the damping tail of the anisotropy 
power spectrum is found. We also show the correction related to second order 
in the deflection angle power spectrum. The resulting modification related to 
this additional higher order correction is at the level of 20%, at most, though 
it becomes significant as I — > 10^, suggesting that, at small angular scales, 
the perturbative approach may be suspect to describe the lensing effect on 
CMB. Still, the approach considered in certain calculations (e.g., Hu 2000) 
is adequate since CMB anisotropics have no significant power at small an- 
gular scales and higher order terms, which are important at small angular 
scales, have a decreasing significance. For precision cosmology with arcminute 
angular-scale CMB power spectrum measurements, however, it may be use- 
ful to include higher order terms since they are above the cosmic variance 
limit of the power spectrum measurements. In addition to these higher order 
corrections, there are other assumptions that can complicate the interpreta- 
tion of CMB lensing data. For example, the CMB power spectrum contains 
an anisotropy contribution from low redshifts through secondary effects such 
as the integrated Sachs- Wolfe effect (Sachs & Wolfe 1967) and the Sunyaev- 
Zel'dovich effect (Sunyaev & Zel'dovich 1980). Since these contributions are 
also lensed, simply using lensed power spectrum with a narrow source at Tq 
for cosmological parameter interpretation may be problematic. This issue is 
beyond the purpose of the present paper and we will leave it for a future study. 

As we will now discuss, due to the smoothness of spatial fiuctuations generated 
by low redshift diffuse backgrounds at IR wavelengths and at 21 cm rest- 
wavelength by neutral Hydrogen rich IGM prior to reionization, the lensing 
modifications, at arcminute angular scales, are not significant. This limits 
their ability for a lensing reconstruction in the same manner as approached 
for CMB. 



2.2 Diffuse Background Clustering Power Spectrum 

Since gravitational lensing defiections modify the clustering power spectrum 
of anisotropies in the diffuse background, we first consider an estimate of the 
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intrinsic spatial clustering signal. In the introduction, we identified several 
potential backgrounds that may become important for weak lensing studies 
including that due to first stars and galaxies in near-IR wavelengths and the 
21 cm background by the general IGM prior to reionization. Instead of mak- 
ing separate estimates of their clustering signals here, we consider a general 
calculation under the assumption that background sources, either in the form 
of first stars or neutral Hydrogen, trace linear fiuctuations with a bias factor 
that can be calculated based on the halo mass in which these sources are ex- 
pected. In the case of first stars, such halos are likely to have temperatures 
of 200 K, if molecular Hydrogen coohng is allowed, or 10^ K, or more. In the 
case of diffuse emission related to, say 21 cm, from the neutral hydrogen prior 
to reionization, we expect neutral Hydrogen to be present in all halos with no 
cut off at the low end of the mass distribution. In this case, it is likely that 
the neutral Hydrogen distribution is tracing the linear density field directly 
such that the bias factor is close to unity. 

The clustering of any sources at redshifts between 10 and 30 can be calculated 
following standard approaches (Cooray et al. 2003). The contribution to the 
background intensity, say at a certain wavelength and towards a direction n, 
can be written as a product of the mean emissivity and its fiuctuation 



/(n) = j dra{r)j{r) 



Sj{rh,r) 



(18) 



where r is the comoving radial distance (or conformal time) and j{r) is the 
mean emissivity per comoving unit volume as a function of distance. Note that 
the integrand can now be identified with 61{rh,r) in Eq. 1. In order to cal- 
culate spatial fiuctuations related to the emissivity, we assume 5j{rh,r)/ j{r) 
trace fiuctuations in the source density field, 5s — 5ps/ps such that, in Fourier 
space, 

^ = «k,.). (19) 



The density field fiuctuations of emitters are defined by the three dimensional 
power spectrum, which we define as 

(5,(k, z)5s{\^, z)) = {2T:f5D{\^ + k')P..(A;, z) . (20) 



To calculate the source power spectrum, as discussed, we scale the linear 
clustering power spectrum by a bias factor. While in the case of 21 cm emission 
from neutral gas the source bias is expected to mostly unity, for sources that 
form in halos above a certain mass cut-off, Mcut , the bias factor is simply given 
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by 



/Mcut dm m n{m) 



(21) 



where 6(m) is the halo bias (Mo & White 1996; Mo et aL 1997) with respect 
to the density field and n{m) is the mass function (Press & Schechter 1974). 

Using the Limber approximation (Limber 1954), the angular power spectrum 
for a distribution of sources that trace a three-dimensional power spectrum 
Pss{k), when projected on the sky, is given by 

Ci = l dT^f{z)G\T){hM{r)fP''- (k = ^,r = oY (22) 



and the integrand is Cf^ir) in Eq. 10. 

The power spectrum of diffuse emission related to sources in IR is illustrated 
in Fig. 3. The calculation related to j{z) in this case follows that of Cooray et 
al. (2003) and we make use of the most optimistic model there for illustrative 
purposes here. The shape of the angular power spectrum is essentially a reflec- 
tion of the projected linear power spectrum over the redshift distribution of 
background sources. Since first stars are more likely to be found in the similar 
redshift range as of the neutral Hydrogen, which contributes to the 21 cm 
background, we expect the shape of the anisotropy clustering power spectrum 
of the 21 cm background to remain the same. Due to differences in emissivity, 
of course, the clustering power spectrum will have a different amplitude in this 
case. 

In Fig. 3, we also show the angular power spectrum of spatial fluctuations in 
the background diffuse emission after gravitational Icnsing. While the lensing 
effect results in a transfer of power from arcminute angular scales, where the 
angular power spectrum peaks, to tens of degree angular scales, where the 
power spectrum is rapidly rising, there is no substantial modification to the 
angular power spectrum at arcminute angular scales. In comparison, as shown 
in Fig. 2, the lensing effect on CMB anisotropics results in a transfer of power 
from degree angular scales, where the power spectrum peaks in this case, to 
arcminute angular scales such that there is a substantial enhancement of power 
along the damping tail by more than two orders of magnitude. The lack of 
lensing modification to the clustering power spectrum of diffuse backgrounds 
can be understood based on Eq. 15 and results from that fact that clustering 
at arcminute scales is smooth and results from the collective emission from a 
large number of sources such that under gravitational lensing, though there 
may be modifications to the distribution of these sources at arcsecond angular 
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scales, the distribution at arcminute angular scales remain the same. If any 
differences exist, these are at the level substantially below a percent. The lack 
of a substantial difference in the angular power spectrum before and after 
gravitational lensing suggests that statistics that primarily target two-point 
clustering information will not be useful for lensing studies. For example, the 
variance related to the distribution of these background fluctuations will not 
be substantially varying and we suggest that the variance measurements may 
not be ideal. The lack of a major difference also limits the extent to which 
lensing information can be extracted from non-Gaussian statistics, such as 
those related to the trispectrum that have been considered with respect lensing 
studies with CMB. 

As discussed earlier, our calculation on the lensing effect on background clus- 
tering power spectrum may be incomplete. Since the typical deflection angle 
is of order an arcminute and there is substantial clustering power at arcminute 
scales, it is necessary that one properly accounts for higher order corrections 
to the previous calculations. We have attempted this numerically and show 
the first and second order contributions in Fig. 3. Unlike the case with CMB 
anisotropics, where the second order contribution was lower than the first or- 
der term, at angular scales of arcminutes and below, we find that this to be 
no longer the case here. In the multipole range of interest, however, we expect 
there to be no substantial contribution from the second and other higher or- 
der terms, though if we were to study the lensing corrections to the clustering 
power spectrum at arcsecond angular scales, then it may be necessary to per- 
form an improved calculation. We restrict to arcminute scale clustering here 
since such clustering may easily be detectable in near future. On the other 
hand, even if clustering at arcsecond angular scales is detected, we expect at 
such angular scales, the shot-noise associated with finite density of sources 
that contribute to the background emission to be the dominant contribution 
to clustering. In Fig. 3, for illustration, we show the expected shot-noise re- 
lated to the IR background generated by first stars and galaxies following 
calculations in Cooray et al (2003). 

While we have specifically discussed the lensing modification to the IR back- 
ground generated by first stars, our conclusions apply to clustering aspects 
of all other backgrounds. As argued earlier, for backgrounds at the same red- 
shift range, one expects a similar shape to the angular power spectrum though 
the amplitude is different for obvious reasons. Our calculation related to the 
lensing modification of the background fluctuation power spectrum is inde- 
pendent of its amplitude and the difference between the lensed and unlensed 
power spectra scales accordingly such that the fractional difference between 
the two remain the same. The only differences are, however, associated with 
secondary considerations such as the extent to which shot-noise related to the 
finite density of sources that contribute to each of these backgrounds may 
become a source of confusion noise for lensing reconstruction. 
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Fig. 2. The lensed CMB temperature anisotropy power spectrum. The curves labeled 
'R' and 'R^' show the lensing calculation up to the first- and second-order in Cf ^ 
respectively. Note that typical calculations in the literature involve the first order 
calculation, while the next-order correction considered here, the R^ contribution, is 
at the level of 30% and below and can be ignored for most practical purposes; for 
precision cosmology, where lensing information is used to extract cosmological pa- 
rameters, such higher order corrections must be accounted for and we illustrate this 
why by showing the cosmic variance limit (line labled 'Cosmic Var.') of the power 
spectrum measurements at each of these multipoles with a full-sky experiment. 

2.3 Lensing Reconstruction from the Diffuse Emission 

As a first approach, one can consider a lensing reconstruction using diffuse 
background anisotropy data similar to that suggested for CMB. The idea 
behind here is the presence of a significant non-Gaussian signal in the form 
of a trispectrum, and an approach to probe this non-Gaussianity based on 
quadratic statistics such the power spectrum of CMB-gradients (Zaidarriaga 
& Seljak 1999; Hu 2001; Hirata & Seljak 2003). Here, we make use of the 
discussion in Cooray & Kesden (2003; see, also, Hu 2001; Hu & Okamoto 
2002). Under gravitational lensing, the trispectrum for diffuse backgrounds 
takes the form of 



T'(li, 12,13,14) 



/ 



^/</Cf/[qtf^,3|(r)(li + I3) • l3(li + I3) ■ I4 




(23) 
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Fig. 3. The clustering power spectrum of a diffuse background at near-IR wave- 
lengths due to first stars and galaxies and at redshifts between 10 and 30. The 
dotted hue shows the intrinsic angular power spectrum, while the correction related 
to gravitational lensing is shown with a dashed line. The solid line is the angular 
power spectrum one observes today; at large angular scales, this power spectrum 
is few percent larger than the intrinsic power spectrum while the differences are 
insignificant at arcminute angular scales. The dashed line shows an estimate on the 
shot-noise contribution to the clustering signal when integrated over the finite num- 
ber density of background sources; for diffuse emissions such as 21cm line from the 
general IGM prior to reionization, we expect this shot-noise to be significantly lower 
or non-existent. The curves labeled 'R' and 'R^' show the first and second-order 
contributions related to the lensing calculation. The curve labeled 'magnification' 
is an estimate on the magnification-related correction to the power spectrum when 
a finite number density of sources is involved and the background is studied after 
removing certain resolved bright sources. 

where the permutations contain 5 additional terms with the replacement of 
(^3; h) pair by other combination of pairs. The derivation related to Eq. 23 can 
be found in Cooray & Kesden (2003), though, we have modified it slightly to 
account for the finite width in the background source distribution by including 
an explicit integral over the radial distribution. The lensing extraction make 
use of this trispectrum based on the fact that it is simply proportional to Cf"^. 

In practice, one makes an estimator for the deflection angle based on quadratic 
statistics. In the present case, when compared to CMB, note that the one esti- 
mates a weighted average of the deflection angle over the redshift distribution 
of sources that contribute to the diffuse emission. The estimator is such that. 
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Fig. 4. The extraction of lensing information from clustering data related to the dif- 
fuse background. The solid lines show the power spectrum of projected potentials 
while the dotted line show the expected noise associated with the lensing recon- 
struction. Here, we show the case at 2; ~ 15, based on, for example, IR background 
from first stars and galaxies, or 21 cm emission from neutral Hydrogen, while, for 
comparison, we also show the case with CMB using, for example, e arcminute scale 
resolution data expected from a mission such as CMBpol. In comparison, CMB 
temperature fluctuations allow a significantly better reconstruction related to grav- 
itational lensing than low redshift backgrounds and is primarily due to the difference 
in the fluctuation pattern on the sky or the clustering power spectrum. The presence 
of features and the lack of power at arcminute scales in CMB anisotropics enhance 
the lensing modification, while in the case of diffuse backgrounds, the lensing mod- 
ifications are minor and the generated non-Gaussianities are no significant. 

in Fourier space, one is convolving two filtered versions of the background 
intensity fluctuations 



where is a filter that is designed to optimally extract information related 
to the weak lensing effect such that the power spectrum of returns (Cf ) 
where 




(24) 




(25) 
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where n.s(r) is the normahzed distribution of background sources that con- 
tribute to the diffuse emission. The fiher, li), can be represented by 
typical quadratic combinations considered in the literature, such as gradients 
(e.g., Zaldarriaga & Seljak 1999). To extract all information, following the case 
in CMB anisotropies, the optimized filter to extract lensing information takes 
the form of 



[i-iiC;, + i-(i-ii)qi_m] 

^'^ii '^ll-lil 



(26) 



such that when one constructs N^{P{\)P{V)) one measures (Cf) with a noise 
contribution given by 



N, 



d\ [1 • hCi, + 1 ■ (1 - ii)qi-u|]^ 



O/^tot/^tOt 



(27) 



Here C*°* is the total contribution to the angular power spectrum of the diffuse 
background map and includes contributions such as C*°* = C/ + Cf + C™^^ + 
Qnoise^ which are related to the unlensed clustering power spectrum including 
any shot-noise, Ci + Cf^°^, any corrections related to the magnification, if again 
there is a finite number density of sources involved C™^^ and instrumental 
noise, Cf™^^ 

The shot-noise associated with a finite density of sources that determine the 
background emission can be estimated through number counts, dN/ dS, of the 
contributing sources, as a function of fiux S, and can be written as 



^shot 



Scut , , , 



(28) 



where S'cut is the flux cut off value related to the removal of resolved sources. 
This shot-noise acts as a source of noise for lensing reconstruction. Unlike the 
calculation related to the clustering spectrum, where the average emissivity 
is the only unknown, the shot-noise depends strongly on detailed aspects of 
source number counts. Since one weighs by a factor S^, the shot- noise is more 
sensitive to rare objects which are brighter. Since there are many uncertainties 
related to establishing the source number counts precisely, we loosely estimate 
this following Cooray et al. (2003) and show this estimate in Fig. 3 for the case 
of IR-background generated by first stars and galaxies. Note that for "truly" 
diffuse emissions, such as the 21 cm rest-wavelength background from general 
IGM and CMB, this is not a concern. 

Note that, so far, we have considered the case where the clustering analysis 
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takes in to account all contributions to the background. In some situations, 
such as in cases where some part of the background is due to resolved sources 
in data, one can study clustering properties by removing such sources from 
the analysis. While we do not expect this to be the case either with near- 
IR background from first objects or in 21 cm emission, most dusty galaxies 
that contribute to the far-IR background can be resolved and removed. In 
such the clustering power spectrum of the residual emission contain a 

correction associated with lensing magnification (Tegmark & Villumsen 1997) 
and can be written as 

(29) 

'S'cut 

This noise contribution depends on the number counts at the flux-cut off. 
If this contribution can be identified, it alone may be used for a lensing 
study though the modification to the power spectrum is expected to be small. 
Though we make an estimate of the contribution in Fig. 3 (curve labeled 
'magnification'), this estimate should also be considered as highly uncertain 
since number counts of sources that contribute to the background are not well 
defined in the context of backgrounds we have primarily considered. This ef- 
fect is again only present in the case of a finite density of sources and is not 
a concern when, for example, all sources that contribute to the background 
have the same flux. 

In addition to effects related to a flnite density of sources that lead to a 
diffuse background, the presence of instrumental noise also affect the lensing 
reconstruction. This noise is given by 

Cr = UyW-'e^'-'^ , (30) 

when /sky is the fraction of sky surveyed, = Airap^^/ Npi^ is the variance 
per unit are on the sky with an individual pixel noise variance of with 
A^pix pixels, and is the effective beam-width of the instrument. 

In Fig. 4, we illustrate the extent to which CMB and low redshift diffuse 
backgrounds can be used for a lensing reconstruction by concentrating on the 
expected level of noise associated with the reconstructed potential power spec- 
trum. For simplicity, here we consider arcminute scale resolution experiments, 
but set the noise level to be ten orders of magnitude below the clustering power 
spectrum such that at multipoles of few hundred, the noise-level is insignif- 
icant. This is done to consider the full extent to which diffuse backgrounds 
can be used for lensing studies based on clustering information and to com- 
pare with CMB related lensing studies. As shown in Fig. 4, CMB data allow 
an order of magnitude or better reconstruction of the projected lensing po- 
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tential power spectrum when compared to the case with diffuse baclcgrounds. 
This is due reasons we have aheady mentioned and involve the fact that the 
non-Gaussianity generated in CMB data, by the nature of its fluctuation pat- 
tern, is more significant than the non-Gaussianity generated by gravitational 
lensing in the low redshift diffuse emission. While fluctuations in the diffuse 
emission may allow a reconstruction of the lensing power spectrum with a 
cumulative signal-to-noise ratio of order few hundred, this should be com- 
pared to at least an order of magnitude better reconstruction with CMB with 
the same angular resolution and fractionally the same noise level. Unlike low 
redshift backgrounds, CMB also has the advantage that a lensing reconstruc- 
tion can be considered from polarization information in addition to tempera- 
ture anisotropics we have only considered here (e.g., Guzik et al. 2001; Hu & 
Okamoto 2002). 

In addition to problems related to the reconstruction, there are also reasons 
that can complicate a simple interpretation of the reconstructed lensing power 
spectrum. For example, in addition to the lensing contributions we have dis- 
cussed, the lensed clustering power spectrum also contains additional contri- 
butions related to the fact that the source distribution that contribute to low 
redshift background is broad and these lensing potentials can be correlated 
with the diffuse emission itself. We write the dominant contribution associ- 
ated with such a cross-correlated component as 

^ / (£y^ / ^^i^t''^^^) / dr2Cfrj;iir2m - li) ■ li]^ . (31) 

This term does not exist in the case of CMB since the potentials that deflect 

CMB photons are disjoint in redshift space from last scattering surface where 
most CMB fluctuations arc generated. In the case of diffuse backgrounds, if 
the source distribution is broad, this term can be interpreted simply as the 
contribution due to lensing potentials at the low end of the source distribution 
that is responsible for the lensing of sources at the other end of the distribution. 
Assuming that both potentials and the sources of background light anisotropy 
trace the same power spectrum P{k), with a bias factor b{k,r2) that accounts 
for any departures from this assumption, this cross power spectrum can be 
written as 

Ts OO 

Cf-''(r.) = -/ k''dkP{k) I drW''^{k,r{)3i{kr{) j dr2b{k,r2)b{k,r2)ji{km) 



A proper accounting of effects such as this will be a necessary aspect of any 
study that attempt to use background source clustering for lensing purposes. 

While a reconstruction with CMB is better, note that one reconstructs the 
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potential power spectrum out to a redshift of 1100. Thus, low redshift back- 
grounds are still useful for cosmological studies since they allow a reconstruc- 
tion out to redshifts of 15 to 30 or so and will bridge the gap between low 
redshift lensing studies, based on galaxy shapes, and eventual lensing studies 
with higher resolution CMB data in the near future. In this respect, we sug- 
gest that one should concentrate further on the aspects of lensing related to 
diffuse backgrounds and how they can be exploited. As we have discussed, the 
use of non-Gaussian information for a lensing reconstruction with low redshift 
backgrounds, in the same manner lensing studies with CMB are proposed, is 
not helpful, though an attempt is still encouraged. To extract lensing infor- 
mation, one should concentrate on the distribution of sources, that contribute 
to the background emission, at arcsecond angular scales. One approach would 
be to use the local anisotropy of the correlation function to obtain an esti- 
mate for the shear (van Waerbeke et al. 1997); this comes from the fact while 
clustering properties of the background is expected to isotropic, in addition to 
non-Gaussianities, gravitational lensing also modify the clustering such that, 
locally, an anisotropic correlation function is generated. The quadrupole mo- 
ment of the local correlation function captures this anisotropy and provides an 
estimate of shear just as ellipticity is used to estimate lensing shear. Another 
approach is to use variations in number counts due to lensing magnification 
related effects (Jain 2002), though this requires backgrounds where sources 
counts are easily modified by magnification due to the nature of source count 
distribution as a function of redshift. While we have primarily discussed the 
use of arcminute scale clustering information and the potential presence of 
non-Gaussianities for a lensing reconstruction, in a future paper, we will re- 
turn to the issue of extracting lensing shear from small angular scale clustering 
properties. 



3 Summciry 



The current weak lensing measurements of the large scale structure are mostly 
related to statistical study of background galaxy ellipticities. We consider a 
possibility to extend lensing studies with intrinsically unresolved sources and 
suggest that clustering properties, at arcminute angular scales, of either the 
point sources or the integrated diffuse emission from these point sources can 
be used for a lensing reconstruction. This is analogous to techniques now de- 
veloped to extract lensing information from cosmic microwave background 
(CMB) temperature and polarization fluctuations and relies on the presence 
of a significant non-Gaussian signal in the background generated by the grav- 
itational lensing modifications to the fiuctuation pattern. Contrary to CMB, 
however, the lensing modification to the diffuse background clustering is min- 
imal and results from the lack of distinct features in the fluctuation pattern. 
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In the case of CMB, a distinct signature or pattern exist through the acoustic 
nature and the lack of intrinsic power at arcminute angular scales where the 
damping-tail is found. We have discussed the extent to which projected matter 
power spectrum be reconstructed with images of the diffuse background gen- 
erated by first stars and galaxies in near-IR wavelengths and the background 
related to 21 cm emission by neutral gas in the general intergalactic medium 
prior to reionization and have shown this to be significantly below the level 
of signal-to-noise ratio one can reach with, in relative terms, CMB data. This 
is due to the lack of significant non-Gaussianities in the diffuse emission that 
is generated by gravitational lensing modifications. Though limited, lensing 
studies with diffuse backgrounds will only be the way to extend the recon- 
struction of the integrated matter power spectrum out to redshifts of 15 to 
30 and to bridge the gap between current and upcoming lensing studies with 
sources at redshifts between 1 and 2 and planned weak lensing studies with 
CMB out to the last scattering surface at a redshift of 1100. As ways to im- 
prove, we suggest further studying a potential lensing reconstruction methods 
based on shear estimation from the local anisotropy of the correlation func- 
tion of the sources that contribute to diffuse backgrounds and statistics that 
are optimized to extract information related to modifications associated with 
lensing magnification. 
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